REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  NO.  0704-0188 


Public  Reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188,)  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (  Leave  Blank) 

2.  REPORT  DATE 

9-8-2006 

3.  REPORT  TYPE  AND  DATES  COVERED 

reprint 

4.  TITLE  AND  SUBTITLE 

5.  FUNDING  NUMBERS 

Effects  of  Alloying  on  the  Chemistry  of  CO  and  Ha  S  on  Fe  Surfaces 

DAA1 9-03-1 -0022 

6.  AUTHOR(S) 

D.E.  Jiang  and  Emily  A.  Carter 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Princeton  University,  Department  of  Mechanical  &  Aerospace  Engineering, Equad, 
Suite  D404a,  Princeton,  NJ  08544-5263 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

10.  SPONSORING  /  MONITORING 

AGENCY  REPORT  NUMBER 

43763.1 6-CH 

1 1 ,  SUPPLEMENTARY  NOTES : 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 

12  a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

12  b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 

13.  ABSTRACT  (Maximum  200  words) 

14.  SUBJECT  TERMS 

15.  NUMBER  OF  PAGES 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OR  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

ON  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540-01-280-5500 

Standard  Form  298  (Rev.2-89) 

Prescribed  by  ANSI  Std.  239-18 
298-102 

J.  Phys.  Chem.  B  2005,  109,  20469-20478 


20469 


Effects  of  Alloying  on  the  Chemistry  of  CO  and  H2S  on  Fe  Surfaces 

D.  E.  Jiang^ 

Department  of  Chemistry  &  Biochemistry,  Box  951569,  University  of  California, 

Los  Angeles,  California  90095-1569 

Emily  A.  Carter* 

Department  of  Mechanical  and  Aerospace  Engineering  and  Program  in  Applied  and  Computational 
Mathematics,  Princeton  University,  Princeton,  New  Jersey  08544-5263 

Received:  May  20,  2005;  In  Final  Form:  August  12,  2005 


Deleterious  gases  such  as  CO  and  H2S  can  cause  degradation  of  steel  by  reacting  with  the  metal  surface. 
Here  we  consider  whether  alloying  the  steel  surface  might  be  able  to  inhibit  these  damaging  surface  reactions 
by  raising  the  barriers  to  molecular  dissociation.  We  employ  first-principles  density  functional  theory  techniques 
to  investigate  the  elementary  reaction  pathways  and  barriers  for  CO  and  H2S  on  FeAl  and  FesSi  surfaces  and 
compare  them  with  pure  Fe  surfaces  (as  a  model  for  steel).  We  find  that  H2S  dissociates  on  iron  surfaces 
much  more  easily  than  CO  does.  Although  FeAl  surfaces  raise  the  barriers  for  H2S  dissociation,  they 
significantly  lower  the  barriers  for  CO  dissociation.  On  the  other  hand,  Fe3Si  surfaces  raise  the  barriers  for 
CO  dissociation,  but  they  are  as  vulnerable  as  Fe  surfaces  to  H2S  dissociation.  Our  findings  suggest  that 
alloying  iron  with  A1  or  Si  is  unlikely  to  simultaneously  increase  its  resistance  to  the  initial  stages  of  chemical 
degradation  by  CO  and  H2S. 


1.  Introduction 

The  deterioration  of  structural  metals  due  to  corrosion  incurs 
great  cost  to  the  economy.  Investigating  the  gas— metal  chem¬ 
istry  involved  will  help  characterize  the  problem  and  ultimately 
help  design  more  corrosion-resistant  structural  materials.  Among 
the  deleterious  gases  for  structural  metals,  carbon  monoxide 
(CO)  and  hydrogen  sulfide  (H2S)  are  among  the  most  notorious. 
CO  is  the  main  culprit  for  two  high-temperature  steel  corrosion 
phenomena:  carburization  and  metal  dusting.  Both  processes 
start  with  CO  adsorption  and  dissociation  on  steel  surfaces, 
followed  by  C  ingress  into  the  steel  lattice.  Carburization 
happens  at  a  relatively  high  temperature  (>800  °C)  with  the 
formation  of  internal  carbides,  which  are  brittle  and  can  cause 
cracks  in  the  steel  substrate.  Metal  dusting  is  generally  observed 
between  350  and  900  °C  with  a  relatively  high  C  activity,  which 
shifts  the  chemical  equilibrium  toward  graphite  formation. 
Although  metal  dusting  is  preceded  by  carbide  formation,  those 
carbides  are  not  stable  and  decompose  to  fine  metal  particles 
and  graphite  (or  other  forms  of  carbon).^  Formation  of  surface 
layers  of  cementite,  FesC,  with  its  lower  melting  point  (com¬ 
pared  to  steel)  may  also  result  in  melting  of  the  steel  surface. 
Those  metal  particles  or  molten  cementite  can  easily  be  washed 
away  by  the  gas  stream,  causing  erosion  of  steel. 

Gases  containing  H2S  attack  structural  metals  such  as  iron 
aggressively.  Due  to  favorable  bulk  thermodynamics  and  the 
relatively  weak  H— S  bond,  the  barriers  for  H2S  dissociation 
on  transition-metal  surfaces  are  usually  very  small.'^“®  This  leads 
to  fast  sulfur  deposition  on  metal  surfaces  and  subsequent  sulfide 
formation.  Sulfur  can  also  cause  the  embrittlement  of  transition 
metals  such  as  Ni  and  Fe.^  Moreover,  since  H2S  deposits  its 
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hydrogen  atoms  on  transition-metal  surfaces,  those  hydrogen 
atoms  can  also  embrittle  metals  such  as  Fe.^’® 

In  some  special  cases  of  metal  erosion/corrosion,  both  CO 
and  H2S  are  present  in  the  gas  phase,  e.g.,  in  the  propellant 
gases  from  a  gun  firing.'**  In  the  making  of  gun  barrels, 
electroplated  Cr  is  used  to  protect  the  steel  substrate.  However, 
the  electroplating  process  produces  Cr  coatings  that  have 
microcracks.  These  microcracks  grow  and  propagate  to  the 
substrate  metal  during  operation,  allowing  those  deleterious 
gases  to  travel  along  those  cracks,  subsequently  attacking  the 
substrate.  Both  iron  sulfide  and  iron  carbide  are  observed  as 
erosion  products.***  Here  we  explore  the  possibility  of  designing 
an  alternative  thin  film  coating  for  steel  with  various  ferrous 
intermetallic  compounds.  To  survive  harsh  conditions,  the 
intermetallic  compound  should  be  stable  from  room  temperature 
to  a  relatively  high  temperature  (e.g.,  1300  °C),  and  the  alloying 
element  should  not  be  a  known  embrittler  of  steel.  Moreover, 
this  element  should  not  be  so  prone  to  reaction  with  CO  or 
S-containing  molecules.  This  element  optimally  would  form  a 
protective  oxide  scale  when  the  ferrous  intermetallic  compound 
is  exposed  to  air.  From  these  considerations,  we  arrived  at  two 
candidate  compounds  to  investigate:  FeAl  and  Fe3Si. 

Industrial  demand  for  lightweight  structural  materials  for  high- 
temperature  applications  has  driven  the  interest  in  the  research 
and  development  of  intermetallic  compounds."  FeAl-  and  Fe3- 
Si-based  alloys  have  been  shown  to  possess  good  corrosion 
resistance  against  carburization  and  sulfidation.'^”*^  This  re¬ 
sistance  is  thought  to  originate  from  the  thin  oxide  scale  formed 
on  the  alloy  surfaces.  We  consider  here  the  possible  use  of  FeAl 
or  Fe3Si  as  a  thin  film  coating  on  steel.  Then,  an  oxide  scale  of 
alumina  or  silica  may  be  able  to  protect  the  alloy  coating  in 
oxidizing  environments.  However,  in  some  industrial  applica¬ 
tions,  the  gas  environment  is  reducing  in  nature,  and  the  oxide 
scale  may  either  be  destroyed  or  not  form  at  all.  Moreover, 
cracks  in  the  oxide  scale  can  also  lead  to  gas  attack  of  the  base 
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alloy.  In  this  work,  we  will  address  how  the  alloy  surfaces  will 
respond  to  the  presence  of  both  CO  and  H2S  under  such 
conditions. 

Although  molecular  adsorption  and  subsequent  surface  chem¬ 
istry  are  well  studied  for  elemental  metal  surfaces,  chemistry 
on  intermetallic  surfaces  is  less  well  explored.  With  first- 
principles  density  functional  theory  (DPT)  techniques,  we  have 
studied  the  adsorption  and  dissociation  of  CO  and  H2S  on  pure 
Fe  surfaces.  In  the  present  work,  we  compare  those  findings 
with  the  chemistry  of  CO  and  H2S  on  the  Fe-alloy  surfaces 
formed  from  the  Fe  intermetallic  compounds  FeAl  and  FesSi. 
We  first  identify  the  site  preference  for  CO,  HS,  H2S,  and  their 
constituent  atoms  on  low-index  surfaces  of  Fe3Si  and  FeAl.  We 
then  characterize  pathways  for  dissociation  of  these  molecules, 
including  activation  barriers,  and  compare  them  with  those  on 
pure  Fe  surfaces. 

The  paper  is  organized  as  follows.  In  section  2,  we  outline 
the  calculational  details.  Results  for  the  bulk  and  surface 
properties  of  FeAl  and  FesSi  are  presented  in  section  3.1.  We 
then  show  the  results  for  adsorption  in  section  3.2  and 
dissociation  in  section  3.3.  The  surface  chemistry  and  effects 
of  alloying  are  discussed  in  section  3.4.  We  summarize  and 
conclude  in  section  4. 

2.  Calculational  Details 

We  perform  first-principles  calculations  based  on  density 
functional  theory  (DFT).'®’'®  The  Vienna  Ab  Initio  Simulation 
Package  (VASP)^®'^'  is  used  to  solve  the  Kohn— Sham  equations 
with  periodic  boundary  conditions  and  a  plane-wave  basis  set. 
Here  we  employ  Blochl’s  all-electron  projector  augmented  wave 
(PAW)  method,^^  as  implemented  by  Kresse  and  Joubert,^^ 
within  the  frozen  core  approximation.  We  use  the  generalized 
gradient  approximation  (GGA)  of  PBE  for  the  treatment  of 
electron  exchange  and  correlation.^'^  Spin-polarized  calculations 
are  employed  for  Fe  and  FesSi  systems  (bulk  and  surface),  while 
nonspin-polarized  ones  are  employed  for  FeAl  systems,  as 
explained  further  in  section  3.1. 

We  use  a  kinetic  energy  cutoff  of  400  eV  for  all  the 
calculations,  which  converges  the  total  energy  to  ~1  meV/atom 
for  the  primitive  unit  cells  of  bulk  Fe,  FeAl,  and  FesSi.  The 
Monkhorst— Pack  scheme^^  is  used  for  the  k-point  sampling. 
The  first-order  Methfessel— Paxton  method^^  is  used  for  the 
Fermi  surface  smearing,  with  a  width  of  0.1  eV  in  order  to  obtain 
accurate  forces.  An  equilibrium  lattice  constant  of  2.834  A  is 
used  for  ferromagnetic  bcc  Fe,  as  we  obtained  earlier  with  a 
converged  k-mesh  of  15  x  15  x  15.^^  This  lattice  constant 
agrees  well  with  the  experimental  value  of  2.86  A.^*  Bulk  and 
surface  properties  for  FeAl  and  FesSi  are  shown  and  discussed 
in  section  3.1. 

To  model  gaseous  CO,  H2S,  HS,  and  C,  O,  H,  and  S  atoms, 
we  place  a  molecule  or  atom  in  a  10-A  cubic  box.  We  perform 
a  nonspin-polarized  calculation  for  CO  and  H2S,  but  spin- 
polarized  calculations  for  open-shell  C,  O,  HS,  H,  and  S,  where 
the  valence  electron  configuration  used  for  C,  O,  and  S  atoms 
is  approximately  the  ground  state  (spin-polarized  DFT  “wave 
functions”  are  slightly  spin-contaminated).  Good  agreement  with 
experiment  is  achieved  for  molecular  geometries,  vibrational 
frequencies,  and  bond  dissociation  energies.*^  '® 

We  employ  a  five-layer  slab  with  a  12- A  thick  vacuum  layer 
between  periodic  images  to  model  all  the  surfaces  and  put 
adsorbates  on  one  side  of  the  slab;  this  produces  a  small  dipole 
due  to  the  dipole  of  the  molecular  adsorbate  itself  and  to  charge 
transfer  between  the  adsorbate  and  the  metal  surface.  However, 
we  did  not  bother  with  a  dipole  correction  to  the  total  energy 


since  it  is  generally  small  (<0.05  eV/cell).  Only  the  top  two 
layers  of  the  five  substrate  layers  are  allowed  to  relax,  together 
with  the  adsorbate  layer.  The  bottom  three  layers  are  kept  fixed 
in  bulk  positions  to  represent  the  semi-infinite  bulk  crystal 
beneath  the  surface.  Allowing  the  middle  layer  of  the  substrate 
to  relax  only  changes  the  total  energy  of  the  slab  by  ~5  meV. 
When  the  maximum  force  acting  on  each  of  the  relaxed  atoms 
drops  below  0.01  eV/A,  the  structural  relaxation  is  stopped.  In 
this  work,  we  study  the  surface  chemistry  of  CO  and  H2S  in 
the  low  coverage  regime  (0  =  0.125—0.25  ML).  We  use  a 
k-mesh  for  a  surface  cell  that  yields  a  k-point  spacing  of  ~0.03 
(2jr/A)  or  smaller  and  that  converges  the  adsorption  energy  to 
within  0.02  eV/adsorbate. 

The  Climbing  Image  Nudged  Elastic  Band  (CI-NEB) 
method^^’^®  is  used  to  locate  minimum  energy  paths  (MEPs) 
and  transition  states.  An  interpolated  chain  of  configurations 
(images)  between  the  initial  and  final  positions  are  connected 
by  springs  and  relaxed  simultaneously  to  the  MEP.  With  the 
climbing  image  scheme,  the  highest-energy  image  climbs  uphill 
to  the  saddle  point.  We  use  the  same  force  tolerance  for  the 
transition-state  search  as  used  for  structural  relaxations.  Typi¬ 
cally,  we  first  use  the  normal  NEB  method  for  ~10  iterations 
(ionic  steps)  to  roughly  converge  the  MEP;  then  we  switch  on 
the  climbing-image  algorithm.  After  ~50  more  iterations,  we 
can  converge  both  the  MEP  and  TS  to  the  preset  force  tolerance. 

Vibrational  frequencies  of  adsorbates  on  surfaces  are  deter¬ 
mined  by  diagonalizing  a  finite  difference  construction  of  the 
Hessian  matrix  with  displacements  of  0.02  A  (only  allowing 
adsorbate  atoms  to  move).  The  natures  of  the  relaxed  adsorbate 
configurations  and  the  saddle  points  found  by  the  CI-NEB 
method  are  characterized  by  vibrational  frequency  analyses, 
where  the  presence  of  one  imaginary  frequency  indicates  a 
transition  state  has  been  found,  while  >2  imaginary  frequencies 
indicates  a  higher-order  saddle  point.  The  absence  of  imaginary 
frequencies  is  associated  with  a  true  minimum  on  the  potential 
energy  surface. 

We  also  use  calculated  frequencies  at  the  adsorption  minima 
and  transition  states  to  estimate  the  zero-point  energies  (ZPE) 
at  those  configurations.  We  found  that  including  ZPE  lowers 
reaction  barriers  by  0.08  eV  or  less.  Similarly,  the  resulting 
changes  to  the  adsorption  energies  due  to  ZPE  are  found  to  be 
<+0.08  eV  (per  CO  or  H2S).  Since  these  changes  are  relatively 
small  compared  to  other  errors  in  the  calculations  and  do  not 
affect  the  general  conclusions  from  the  present  work,  we  do 
not  report  ZPE  corrections  to  the  energetics. 

3.  Results  and  Discussion 

3.1.  Bulk  and  Surface  Properties  of  FeAl  and  FesSi.  The 

ground  state  of  bulk  Fe  is  the  body-centered-cubic  (bcc) 
ferromagnetic  (FM)  phase.  In  our  earlier  work,  we  explored  the 
various  phases  of  Fe  with  PAW-DFT-GGA  and  compared  the 
results  with  experiment.^^  Good  agreement  is  obtained  for  the 
lattice  parameter,  local  magnetic  moment,  and  bulk  modulus 
of  bcc  FM  Fe.  The  structures  of  FeAl  and  FesSi  are  related  to 
bcc  Fe.  We  present  their  structural  properties  next. 

3.1.1.  Bulk  FeAl  and  FejSi.  FeAl  has  a  CsCl  (B2)  type  of 
structure,  that  is,  one  A1  atom  replacing  Fe  at  the  cubic  center 
of  bcc  Fe  (see  Figure  1).  Although  magnetic  susceptibility, 
Mossbauer,  and  nuclear  magnetic  resonance  measurements 
indicate  that  perfectly  ordered,  stoichiometric  FeAl  is  nonmag¬ 
netic  in  nature,^ an  FM  ground  state  is  found  by  DFT  (both 
LDA  and  GGA)  with  a  local  magnetic  moment  of  ~0.7  for 
Fe.^^  This  FM  state  is  only  slightly  lower  in  energy  (7  meV/ 
atom)  than  the  nonmagnetic  (NM)  state.  Since  the  FM  and  NM 


Chemistry  of  CO  and  H2S  on  Fe  Alloy  Surfaces 


J.  Phys.  Chem.  B,  Vol.  109,  No.  43,  2005  20471 


Figure  1.  Bulk  and  surface  structures  of  Fe,  FeAl,  and  FesSi.  Fe  atoms 
ai'e  in  black,  A1  in  gray,  and  Si  in  white. 


TABLE  1:  Comparison  of  DFT  Predictions  with 
Experiment  for  Bulk  Properties  of  FesSl:  Lattice  Parameter 
(a),  Bulk  Modulus  (B),  and  Local  Magnetic  Moments  (M) 

Mpe-B  Mpe-A  Msi 

a  (A)  B  (GPa)  (mb)  (ub)  (/<b) 

PAW-DFT-GGA-PBE  5.61  194  2.53  1.31  -0.05 

experiment  5.65“  187''  2.2-2.4“  1.35“  -0.07“ 

“Ref  38.  ''Ref  39.  “Ref  37. 

states  are  almost  degenerate  and  there  are  still  controversies 
regarding  the  true  ground  state  of  FeAl,^'*’^^  we  choose  the  NM 
state  to  model  bulk  FeAl  and  surfaces  for  simplicity.  We  obtain 
an  FeAl  equilibrium  lattice  constant  of  2.87  A,  which  agrees 
well  with  previous  all-electron  FLAPW-DFT-LDA  calculations 
(2.87  A)  and  experiment  (2.91  A),  and  an  FeAl  bulk  modulus 
of  172  GPa,  which  improves  upon  the  FLAPW-DFT-LDA  result 
(200  GPa)  compared  with  experiment  (152  GPa).^®  The 
improved  bulk  modulus  is  no  doubt  due  to  the  use  of  the  GGA 
for  exchange  and  correlation.  These  results  for  bulk  FeAl  suggest 
that  PAW-DFT-GGA  should  provide  a  fair  description  of  FeAl 
surfaces. 

Bulk  FesSi  has  a  DOs-type  of  structure.  It  can  be  visualized 
as  a  (2  X  2  X  2)  supercell  of  bcc  Fe  with  four  (out  of  eight) 
body-center  Fe  atoms  replaced  by  Si.  Those  four  Si  atoms  form 
a  tetrahedron  (Figure  1).  We  see  that  there  are  two  types  (A 
and  B)  of  Fe  atoms  in  the  bulk  FeaSi  structure:  the  A-type  is 
surrounded  by  four  Fe  atoms  and  four  Si  atoms  in  the  first 
coordination  shell,  while  the  B-type  is  surrounded  by  eight 
nearest-neighbor  Fe  atoms.  Magnetic  susceptibility  measure¬ 
ments  show  that  FesSi  is  FM.^^  We  compare  our  results  for 
bulk  FesSi  with  experiment^^“^^  in  Table  1  for  the  lattice 
parameter,  bulk  modulus,  and  local  magnetic  moments.  The 
agreement  is  quite  good,  suggesting  that  since  PAW-DFT-GGA 
describes  bulk  FesSi  well,  the  surfaces  of  FesSi  may  also  be 
reasonably  well-described  by  this  level  of  theory. 


Figure  2.  High- symmetry  adsorption  sites  on  (a)  Fe(llO):  1  =  OT, 

2  =  SB,  3  =  TF,  and  4  =  LB;  (b)  FeAl(l  10):  1  =  OT-Fe,  2  =  OT-Al, 

3  =  SB,  4  =  TF-2Fe,  5  =  TF-2A1,  6  =  LB-2Fe,  and  7  =  LB-2A1;  (c) 
FejSi)!  10):  1  =  OT-FeA,  2  =  OT-Si,  3  =  OT-FeB,  4  =  SB-2Fe,  5  = 
SB-FeSi,  6  =  TF-FeSi,  7  =  TF-2Fe,  8  =  TF-3Fe,  9  =  LB-4Fe,  10 
=  LB-2Fe,  and  11  =  LB-FeSi.  (Large  circles  are  surface  atoms,  and 
small  ones  indicate  adsorption  sites.  For  Fe3Si,  the  A-type  Fe  is  denoted 
as  a  white  circle  and  the  B-type  as  a  black  circle.  OT  =  on-top,  SB  = 
short-bridge,  LB  =  long-bridge,  and  TF  =  3-fold.) 

3.1.2.  FeAl  and  FesSi  Surfaces.  We  consider  low  Miller-index 
surfaces,  since  they  are  the  most  likely  to  be  exposed.  Because 
FeAl  and  FesSi  both  belong  to  the  bcc  framework,  (110)  will 
be  the  closest-packed  surface.  The  surface  energy  of  Fe(lOO)  is 
predicted  by  DFT-GGA  to  be  only  slightly  higher  than  Fe- 
(1 10),'^°  so  we  also  consider  (100)  surfaces  for  FeAl  and  FesSi. 
Figure  1  displays  bulk- terminated  (110)  and  (100)  surfaces  for 
Fe,  FeAl,  and  FesSi.  Although  recent  surface  science  studies 
(with  low-energy  electron  diffraction  and  Auger  electron 
spectroscopy)  show  that  A1  and  Si  tend  to  segregate  on  Fe— 
intermetallic  surfaces  after  annealing,^'"'**  our  work  focused  on 
the  bulk-terminated  surfaces  for  simplicity.  We  note  that  some 
bulk-terminated  surfaces  are  indeed  obtained  experimentally 
under  certain  conditions  for  FeAl  and  FesSi.  For  example,  bulk- 
terminated  FeAl(l  10)  can  be  obtained  by  annealing  the  surface 
at  400  while  bulk-terminated  FeAl(lOO)  with  A1  at  the 
surface  can  be  obtained  by  annealing  at  927  We  therefore 
use  the  Al-terminated  FeAl(lOO)  in  our  work.  This  surface  is 
expected  to  have  a  lower  surface  energy  than  the  Fe-terminated 
structure,  consistent  with  the  lower  surface  energy  of  pure  A1 
compared  to  pure  Fe.  For  FesSi(lOO),  in  principle  the  surface 
could  be  terminated  with  either  100%  Fe  atoms  or  50%  Fe  and 
50%  Si.  Because  experiment  shows  that  Si  tends  to  segregate 
to  the  surface,"*^  we  investigate  here  the  50:50  terminated  surface 
of  FesSi(lOO),  as  shown  in  Figure  1.  FeAl(llO)  and  FesSi(l  10) 
each  have  only  one  bulk-terminated  structure,  which  are  also 
shown  in  Figure  1. 

With  DFT-GGA-PBE,  we  explored  the  structural  relaxation 
of  bulk-terminated  FeAl  and  FeaSi  surfaces.  We  find  that  the 
spacing  between  the  surface  A1  layer  and  subsurface  Fe  layer 
of  FeAl(lOO)  contracts  by  13.7%,  in  good  agreement  with 
experiment  (14.6%).^°  We  also  find  a  slight  buckling  (~0.1  A) 
in  the  surface  layers  of  FeAl(llO),  Fe3Si(110),  and  Fe3Si(100). 

3.2.  Adsorption  Geometries  and  Energetics.  High-sym- 
metry  adsorption  sites  on  the  (110)  and  (100)  alloy  surfaces 
are  shown  in  Figures  2  and  3,  respectively.  From  Fe(llO)  to 
FeAl(l  10)  to  Fe3Si(l  10),  the  number  of  distinct  sites  increases; 
we  identified  11  sites  on  Fe3Si(110)  to  explore. 

3.2.1.  Adsorption  on  (110)  Surfaces.  Table  2  displays  the  most 
stable  sites  and  the  adsorption  energies  for  all  adsorbates 
considered  on  all  three  (110)  surfaces.  We  predict  that  all  open- 
shell  species  bind  strongly  (due  to  covalent  bonding),  while 
closed-shell  species  are  less  strongly  bound  (due  to  weaker 
donor/acceptor  bonding)  to  these  surfaces.  We  see  that  the  C 
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Figure  3.  High-symmetry  adsorption  sites  on  Fe(lOO),  FeAl(lOO),  and 
Fe3Si(100):  4-fold  hollow  (H),  2-fold  bridge  (B),  and  1-fold  on-top 
(T)  site.  For  Fe3Si(100),  the  on-top  sites  are  divided  into  two  types: 
1-fold-Fe  and  1-fold-Si. 


TABLE  2:  Most  Stable  Adsorption  Sites  and  Energies  (in 
eV t  Ead  EA/Metal— slab  Ejvietal— slab  Ea)  nf  C,  O,  CO,  H,  S, 

HS,  and  HjS  on  Ee(llO),  FeAl(llO),  and  FeaSiillO)" 


species 

properties 

Fe(llO) 

FeAl(llO) 

FeaSilllO) 

c 

site 

LB 

LB-2Fe 

LB-4Fe 

End 

-7.8 

-8.2 

-7.7 

0 

site 

TF 

TF-2Fe 

LB-4Fe 

Eiid 

-6.3 

-6.5 

-6.2 

CO 

site 

OT 

OT-Fe 

LB-4Fe 

End 

-1.9 

-2.5 

-2.2 

H 

site 

TF 

TF-2Fe 

TF-3Fe 

End 

-3.0 

-3.1 

-3.2 

s 

site 

LB 

LB-2Fe 

LB-4Fe 

End 

-5.8 

-5.8 

-5.8 

HS 

site 

LB 

LB-2Fe 

LB-4Fe 

End 

-3.5 

-3.8 

-3.6 

HjS 

site 

SB 

LB-2Fe 

LB-4Fe 

End 

-0.6 

-1.2 

-0.9 

“  Site  definitions  are  given  in  Figure  2. 


Figure  4.  Tilted  states  of  CO  on  alloy  surfaces:  (a)  FeAl(llO);  (b) 
Fe3Si(l  10);  (c)  FeAl(lOO).  Fe  atoms  in  dark  gray,  C  in  black,  and  O  in 
white,  and  A1  and  Si  atoms  in  light  gray. 

atom  binds  strongly  with  all  three  surfaces  and  prefers  high- 
coordination  sites  that  maximize  the  number  of  Fe  atoms.  The 
site  preference  of  S  and  HS  on  (110)  surfaces  is  the  same  as 
that  of  C,  only  the  binding  is  less  strong  than  for  C.  O  also 
binds  strongly  to  (110)  surfaces.  As  on  Fe(llO),  O  prefers  the 
3-fold  (TF)  site  on  FeAl(llO),  but  it  prefers  the  long-bridge- 
4Fe  (LB-4Fe)  site  on  FesSifllO).  Indeed,  the  TF-3Fe  site  on 
Fe3Si(l  10)  is  not  stable;  the  O  atom  relaxes  to  the  LB-4Fe  site 
when  it  is  initially  put  at  the  TF-3Fe  site.  The  H  atom  prefers 
the  TF  site  that  provides  the  maximum  number  of  Fe  atoms, 
and  it  binds  strongly  with  all  three  (110)  surfaces. 

CO  prefers  the  on-top  (OT)  site  on  Fe(llO)  in  an  upright 
fashion.  This  is  also  the  case  for  CO  on  FeAl(l  10),  where  CO 
strongly  adsorbs  in  a  linear  fashion  on  top  of  an  Fe  atom,  while 
binding  of  CO  to  an  A1  site  on  FeAl(llO)  is  ~2  eV  weaker. 
We  also  find  a  tilted  state  of  CO  on  FeAl(l  10)  that  is  ~0.2  eV 
less  stable  than  the  upright  OT  configuration.  We  show  this 
structure  in  Figure  4a,  where  we  see  that  the  C  atom  is  displaced 
by  ~1.0  A  off  the  OT  site  along  [001],  and  the  molecular  axis 
tilts  off  the  surface  normal  (toward  [001]  by  35°).  The  CO 
stretching  frequency  is  determined  to  be  1762  cm“'.  Interest- 


TABLE  3:  Most  Stable  Adsorption  Sites  and  Energies  (in 
eV ;  Ead  E^y^xatal— slab  E^xgtal— slab  E4)  of  C,  O,  CO,  H,  S, 

HS,  and  HjS  on  Ee(lOO),  EeAI(lOO),  and  Ee3Si(100)'' 


species 

properties 

Fe(lOO) 

FeAl(lOO) 

Fe3Si(100) 

c 

site 

4-fold 

4-fold 

4-fold 

End 

-8.2 

-6.8 

-6.7 

0 

site 

4-fold 

2-fold 

2-fold 

Eld 

-6.4 

-6.5 

-6.1 

CO 

site 

4-fold-tilt 

1-fold 

1-fold-Fe 

Ead 

-2.1 

-0.8 

-1.2 

H 

site 

4-fold 

1-fold-Al 

1-fold-Si 

Ead 

-2.7 

-2.4 

-2.6 

S 

site 

4-fold 

4-fold 

2-fold 

Ead 

-6.0 

-4.9 

-4.5 

HS 

site 

4-fold 

2-fold 

2-fold 

Ead 

-3.6 

-3.3 

-2.9 

H,S 

site 

2-fold 

1-fold 

1-fold-Fe 

Ead 

-0.5 

-0.2 

-0.5 

“  Site  definitions  are 

given  in  Figure  3. 

ingly,  a  similar  tilted  state  of  CO  is  also  found 

on  Fe3Si(l  10) 

(Figure  4b).  Here  the  displacement  off  the  OT  site  is  0.73  A, 

the  tilt  angle  is  25°,  and  the  CO  stretching  frequency  is  1900 
cm“*.  Although  CO  prefers  to  reside  in  the  LB-4Fe  site  in  an 
upright  fashion  on  Fe3Si(l  10),  the  tilted  state  is  almost  degener¬ 
ate  with  the  LB-4Fe  site  (within  0.02  eV/CO). 

The  interaction  between  H2S  and  metal  surfaces  is  weakest 
among  the  adsorbates  explored  here,  because  H2S  does  not  have 
either  unpaired  electrons  or  low-lying  empty  jr-orbitals.  It  can 
only  interact  with  the  surface  via  its  lone  pairs  in  a  cr-donor 
mode.  We  see  that  H2S  also  prefers  binding  to  pure  Fe  sites, 
and  it  binds  more  strongly  to  FeAl(l  10)  and  Fe3Si(l  10)  than  to 
Fe(llO).  Electron  density  changes  upon  adsorption  show  that 
accumulation  of  electron  density  between  the  S  atom  in  H2S 
and  its  two  nearest-neighbor  Fe  atoms  increases  with  the 
concentration  of  alloying  elements  on  (110)  surfaces,  namely, 
going  from  Fe(llO)  to  Fe3Si(110)  to  FeAl(llO).  The  charge 
rearrangement  in  each  case  involves  simply  a  polarization  of 
the  electrons  of  H2S.  This  trend  of  increased  polarization  as 
the  alloying  element  concentration  increases  correlates  well  with 
the  change  in  adsorption  energy.  Thus,  alloying  Fe  appears  to 
increase  the  extent  of  cr-donation  by  H2S. 

3.2.2.  Adsorption  on  (100).  FeAl(lOO)  is  terminated  with  Al, 
while  Fe3Si(100)  is  terminated  with  half  Fe  atoms  and  half  Si 
atoms,  as  shown  in  Figures  1  and  3.  The  Si  atoms  at  the  Fe3- 
Si(lOO)  surface  exhibit  a  near-zero  (—0.07  /tb)  local  magnetic 
moment,  as  they  do  in  the  bulk,  indicating  that  there  are  no 
dangling  bonds  associated  with  them.  Table  3  displays  the  most 
stable  sites  and  their  adsorption  energies  on  (100)  surfaces.  Once 
again,  the  open-shell  species  (atoms  and  SH)  are  much  more 
strongly  bound  than  the  closed-shell  molecules,  due  to  the 
former’s  ability  to  form  strong  covalent  bonds  to  the  metal 
surface.  We  see  that  C  always  prefers  a  4-fold  site  and  adsorbs 
most  strongly  on  Fe(lOO).  The  presence  of  Si  on  the  Fe3Si- 
(100)  surface  destabilizes  C  by  ~1.5  eV;  likewise,  C  adsorption 
on  the  fully  Al- terminated  FeAl(lOO)  surface  is  ~1.4  eV  less 
strongly  bound  than  on  a  pure  Fe(lOO)  surface.  This  is  to  be 
contrasted  with  O  adsorption  on  the  three  surfaces,  where  O 
binds  with  nearly  equal  affinity  to  all  of  them.  This  is  likely 
due  to  the  high  affinity  of  Si  or  Al  for  O.  Moreover,  we  see 
that  O  prefers  the  2-fold  bridge  site  on  FeAl(lOO)  and  Fe3Si- 
(100).  Since  FeAl(lOO)  is  terminated  with  an  Al  layer,  it  makes 
sense  to  compare  it  with  Al(lOO).  Extended  appearance-potential 
fine  structure  analysis^*  shows  that  O  prefers  the  4-fold  hollow 
site  and  sits  slightly  below  the  surface  layer  of  Al(lOO).  Recent 
DFT-embedded  cluster  calculations^^  confirm  this  site  preference 
of  O  atoms  on  Al(lOO).  The  difference  between  Al(lOO)  and 
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(a)  (b) 


Figure  5.  Most  stable  adsorption  sites  for  H  on  (100):  (a)  Fe3Si;  (b) 
FeAl.  Fe  atoms  in  dark  gray,  H  in  white,  and  A1  and  Si  atoms  in  light 
gray. 

FeAl(lOO)  for  O  adsorption  is  likely  due  to  the  fact  that  the 
interlayer  spacing  of  Al(lOO)  is  much  larger  (by  0.6  A)  than 
that  of  FeAl(lOO).  As  a  result,  O  atoms  can  burrow  more  easily 
into  the  surface  layer  at  the  hollow  site,  in  preparation  to  go 
into  the  subsurface  and  to  begin  to  form  aluminum  oxide.  At 
the  2-fold  site  of  Fe3Si(100),  O  interacts  more  with  Si  than  Fe, 
as  evidenced  by  bond  lengths  that  are  shorter  (longer)  than  the 
sum  of  atomic  radii  for  O— Si  (O— Fe):  ro-si  =  1.65  A  and 
ro-Fe  =  2.00  A,  compared  with  ro  -F  rsi  =  1.83  A  and  ro  +  rpe 
=  1.93  A.^^  Oxygen  atoms  tend  to  form  covalent  bonds  with 
Si  atoms,  as  in  the  bulk  Si02,  where  an  O  atom  bridges  two  Si 
atoms.  As  a  result,  it  makes  sense  that  O  atoms  move  to  a  lower 
coordination  2-fold  bridge  site  on  Fe3Si(100),  since  O  atoms 
prefer  to  bond  more  covalently  to  the  Si  atoms  on  the  surface. 

We  find  that  H  prefers  the  4-fold  hollow  site  on  Fe(100),^‘' 
in  agreement  with  experiment.^^  The  2-fold  bridge  site  is  also 
a  minimum  for  FI,  but  the  1-fold  OT  site  is  a  higher-order  saddle 
point.^'*  Unexpectedly,  we  find  that  Ft  prefers  the  1-fold  site  on 
FeAl(lOO)  and  also  the  1-fold  Si  site  of  Fe3Si(100).  This  is  in 
contrast  with  Al(lOO),  where  Ft  prefers  the  2-fold  bridge  site.^® 
Figure  5  shows  the  two  structures.  It  is  rare  for  H  to  adsorb  on 
the  1-fold  site  of  a  metallic  surface  with  4-fold  symmetry. 
This  type  of  bonding  indicates  that  the  H— A1  and  H— Si 
interactions  are  likely  very  covalent  in  character,  resembling, 
for  example,  the  Si— H  bonds  on  the  Si(lOO)  surface,  despite 
the  lack  of  dangling  bond  character  on  the  bare  Fe3Si(100) 
surface.  The  bond  lengths  for  H— A1  and  FI— Si  are  found  to  be 
1.627  and  1.524  A,  respectively,  which  are  similar  to  gaseous 
AlH  (fh-a  1=  1.648  A)  and  SiH4  (fh-si  =  1.480  A),^^  lending 
credence  to  the  idea  that  the  bonding  of  H  to  FeAl(lOO)  and 
Fe3Si(100)  involves  localized,  covalent  bonds.  Site-projected, 


E  -  E,  (eV) 

(a) 


orbital-resolved  local  densities  of  states  (LDOS)  shown  in  Figure 
6  also  support  this  picture.  FI  Is  states  overlap  strongly  with  Si 
3s  states  around  —9  eV  and  3pz  states  between  —6  and  —4  eV 
for  FI/FesSijlOO)  (Figure  6a).  On  FeAl(lOO),  H  Is  states  overlap 
with  A1  3s  between  —7  and  —4  eV  and  3pz  between  —4  and 
—  1  eV  (Figure  6b).  The  corresponding  metal— H  antibonding 
states  appear  above  the  Fermi  level,  a  clear  signature  of  covalent 
bonding  (as  in  the  a/a*  levels  in  H2). 

CO  adsorbs  most  strongly  on  Fe(lOO),  preferring  a  4-fold 
tilted  configuration  over  the  OT  site  by  0.6  eV  at  0.25  ML. 
This  site  was  identified  previously  by  several  different  surface 
science  techniques.  A  very  low  CO  stretching  frequency  (1210 
cm“')  was  observed^®  by  high-resolution  electron  energy  loss 
spectroscopy  (FIREELS).  Our  work  also  yields  a  very  low 
frequency  of  1189  cm“*,  consistent  with  experiment.  For  CO 
adsorption  on  FeAl(lOO),  we  also  find  a  tilted  state  that  is 
essentially  degenerate  (within  0.01  eV)  with  the  OT  site.  At 
both  sites,  CO  is  weakly  bound  to  the  surface.  This  adsorption 
energy  is  in  line  with  the  experimental  observation  that  CO 
physisorbs  on  Al(lOO)  at  low  temperatures  and  desorbs  at  room 
temperature.*  The  tilted  state  of  CO  on  FeAl(lOO)  is  found  to 
have  a  dramatically  lower  CO  stretching  frequency  (853  cm“*) 
than  on  Fe(lOO),  and  the  CO  bond  is  very  much  elongated  (rc-o 
=  1.475  A;  0tiit  =  61.2°;  see  Figure  4c),  suggesting  that  despite 
its  weak  adsorption  energy,  it  is  highly  activated.  This  surprising 
finding  can  be  understood  as  a  balance  between  strong  intrinsic 
adsorption  and  destabilization  of  the  CO  molecule  by  breaking 
its  jT-bonds.®'  On  Fe3Si(100),  CO  prefers  to  adsorb  on  Fe  with 
an  upright,  on-top  fashion.  Although  the  adsorption  site  for  CO 
on  Fe3Si(100)  has  not  been  established  experimentally,  photo¬ 
emission  spectra  imply  that  CO  is  indeed  upright  on  Fe3Si- 
(100).®^  For  comparison,  we  note  HREELS  and  temperature- 
programmed  desorption  (TPD)  results  suggesC^’^'^  that  CO 
adsorbs  on-top  of  an  Si  atom  in  an  upright  fashion  on  Si(lOO), 
with  an  adsorption  energy  of  —0.63  eV  predicted  by  DET  with 
a  hybrid  exchange-correlation  functional  and  a  cluster  model. 
In  our  work,  CO  adsorbs  on  the  Si  atom  of  Fe3Si(100)  with  an 
energy  of  only  —0.4  eV,  much  weaker  than  on  the  Fe  atom, 
indicating  that  Si  repels  CO  on  Fe3Si(100). 

S  adsorbs  most  strongly  on  Fe(lOO),  followed  by  FeAl(lOO), 
and  then  Fe3Si(100).  While  S  prefers  the  high  coordination  sites 
on  the  metallic  surfaces  [Fe(lOO)  and  FeAl(lOO)],  it  switches 
to  the  2-fold  site  on  Fe3Si(100),  as  in  the  case  for  oxygen. 
Generally,  then,  the  presence  of  Si  on  the  Fe  surface  forces  more 
covalent  bonding  to  adsorbates.  FIS  adsorption  follows  the  same 
trend  as  S  in  terms  of  site  preference  and  surface  affinity, 
because  the  adsorption  behavior  of  HS  on  (100)  surfaces  is  also 
dictated  by  the  unpaired  electron  on  S.  The  binding  energies 


Figure  6.  Site-projected,  orbital-resolved  LDOS  for  (a)  the  FI  adatom  and  the  surface  Si  atom  connected  to  the  H  atom  from  Figure  5a;  and  (b) 
the  H  adatom  and  the  surface  A1  atoms  connected  to  the  FI  atom  from  Figure  5b. 
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barriers  (eV) 

Fe(llO) 

Fe(lOO) 

FeAl(llO) 

FeAl(lOO) 

Fe3Si(110) 

FejSiilOO) 

CO  —  C  -f  0 

1.5 

1.1 

1.3 

0.4 

2.4 

1.3 

H2S  —  H  -f  HS 

0.1 

0.3 

0.4 

0.3 

0.1 

0.2 

HS  —  H  -I-  S 

0. 

0.3 

0.3 

0.7 

0. 

0.4 

between  HS  and  (100)  surfaces  are  >1.6  eV  less  strong  than 
between  S  atoms  and  (100)  surfaces,  due  to  the  fewer  number 
of  unpaired  electrons  in  HS.  Interestingly,  HS  is  bonded  linearly 
to  the  Fe(lOO)  surface  at  the  4-fold  hollow  site,  whereas  HS  is 
bent  toward  the  surface  on  both  FeAl(lOO)  and  Fe3Si(100)  at 
the  2-fold  bridge  site  by  66.6°  and  70.0°,  respectively,  indicating 
once  again  that  alloying  the  surface  enhances  covalent  bonding 
to  adsorbed  species. 

H2S  adsorbs  extremely  weakly  on  all  three  surfaces,  with  the 
weakest  interaction  between  H2S  and  FeAl(lOO).  On  both  FeAl- 
(100)  and  Fe3Si(100),  the  H2S  molecular  plane  is  almost  parallel 
to  the  surface.  The  difference  is  that  the  S  atom  of  H2S  binds 
to  A1  on  FeAl(lOO)  but  to  Fe  on  Fe3Si(100).  This  type  of  parallel 
on-top  configuration  is  also  found  for  H2S  on  Fe(lOO),  but  its 
adsorption  energy  is  predicted  to  be  ~0. 1  eV  weaker  than  the 
2-fold  bridge  site,  where  the  molecular  plane  is  oriented 
perpendicular  to  the  surface. 

3.3.  Dissociation  Pathways  and  Transition  States.  After 
the  preferred  adsorption  sites  for  CO,  H2S,  HS,  and  atomic 
species  were  identified  on  the  surfaces,  we  then  located 
minimum  energy  paths  and  transition  states  for  CO  dissociation 
and  H2S  dehydrogenation.  Table  4  displays  the  barriers  we 
obtained  on  the  six  surfaces  for  the  three  dissociation  processes. 

3.3.1.  CO  Dissociation.  After  examining  all  the  structures  at 
the  transition  state  for  CO  dissociation,  we  find  that  they  share 
the  common  feature  that  CO  lies  side-on  at  a  bridge  site, 
presumably  due  to  the  need  to  break  the  jr-bonds  of  CO.  This 
transition-state  configuration  has  been  predicted  previously  for 
CO  and  N2  dissociation  on  metal  surfaces. We  predicted 
that  CO  dissociation  on  Fe(l  10)  has  a  barrier  of  1.5  eV,'®  while 
we  determined  the  barrier  on  Fe(lOO)  to  be  1.1  eV,  which  agrees 
very  well  with  experiment  (1.09  eV)®*  and  previous  DFT 
work.®^’^®  The  lower  barrier  on  Fe(lOO)  can  be  attributed  to  the 
highly  tilted  state  of  CO  on  Fe(lOO),  which  has  a  much 
weakened  CO  bond  already  before  dissociation.  The  OT  site 
preference  by  CO  on  Fe(l  10)  then  requires  more  rearrangement 
of  the  CO  molecule  to  reach  the  TS,  leading  to  a  higher  barrier. 

The  barrier  for  CO  dissociation  on  FeAl  is  found  to  be  lower 
than  that  for  the  corresponding  pure  Fe  surface.  This  is  especially 
true  for  FeAl(lOO).  We  predict  a  barrier  of  only  0.4  eV  on  FeAl- 
(100),  which  is  ~0.7  eV  lower  than  on  pure  Fe(lOO).  Figure  7 
shows  the  minimum  energy  path  (MEP).  This  dissociation 
pathway  is  the  same  as  that  on  Fe(lOO).  The  lower  barrier  on 
FeAl(lOO)  again  can  be  attributed  to  the  much  weakened  initial 
state,  as  we  discussed  above  (Figure  4c).  In  this  state,  the  CO 
bond  length  (1.475  A)  is  even  larger  than  that  in  gaseous 
methanol  (H3C— OH,  1.425  A^*),  indicating  that  the  C— O  bond 
order  in  this  site  is  reduced  to  <  1 .  An  interaction  between  C 
and  the  subsurface  Fe  atom,  as  well  as  the  high  affinity  of  O  to 
surface  A1  atoms,  lead  to  this  highly  tilted  state.  This  is 
supported  by  the  fact  that  we  also  found  this  tilted  state  of  CO 
on  just  a  monolayer  of  A1  on  Fe(lOO),  and  we  demonstrated  a 
significant  interaction  between  subsurface  Fe  and  C  for  this 
latter,  analogous  surface.^'  The  smaller  barriers  on  FeAl  surfaces 
indicate  that  FeAl  alloys  may  be  even  more  susceptible  to 
carburization  than  is  Fe  itself. 

CO  is  found  to  have  a  higher  dissociation  barrier  on  Fe3Si 
than  on  the  corresponding  pure  Fe  surface.  For  example,  the 
CO  dissociation  barrier  is  2.4  eV  on  FesSifllO),  0.9  eV  higher 
than  on  Fe(llO).  The  structure  of  the  MEP  for  breaking  the 


CO  bond  on  Pe3Si(110)  is  found  to  be  similar  to  what  we 
predicted  for  CO  dissociation  on  Pe(llO).'®  Namely,  CO  first 
moves  off  the  OT  site,  then  tilts  toward  the  surface,  and 
dissociates  over  a  bridge  site.  This  pathway  is  also  adopted  by 
CO  dissociation  on  Pe3Si(100),  where  CO  prefers  upright  on- 
top  adsorption  on  the  Fe  atom.  So  CO  dissociation  on  FesSi- 
(100)  is  quite  different  from  that  on  Fe(lOO)  or  FeAl(lOO),  where 
a  side-on  adsorption  state  leads  to  a  smaller  dissociation  barrier. 
The  increase  in  dissociation  barrier  (relative  to  pure  Fe)  on  Fes- 
Si(lOO)  is  smaller  (~0.2  eV)  than  forFe3Si(110).  However,  since 
Fe3Si(110)  has  a  surface  energy  ~10%  less  than  that  of  Fes- 
Si(lOO),  Fe3Si(110)  is  more  likely  to  be  the  preferred  facet 
exposed.  So  overall  Fe3Si  is  expected  to  be  more  resistant  to 
initial  CO  attack  (i.e.,  adsorption  and  dissociation)  than  is  pure 
Fe. 

3.3.2.  H2S  Dissociation.  The  closest-packed  Fe(llO)  surface 
is  found  to  be  very  efficient  in  dissociating  H2S.  The  first 
dehydrogenation  step  has  a  small  barrier  (0.1  eV),  and  the 
second  step  is  barrierless. On  Fe(lOO),  the  barriers  are  slightly 
higher.'^  The  MEPs  show  that  the  H— S  bond  breaks  by  bending 
the  bond  toward  the  surface  and  breaking  off  the  H  atom  on 
top  of  an  Fe  atom.'^’'^  From  Table  4,  we  see  that  the  barriers 
for  the  two  dehydrogenation  steps  on  Fe3Si(110)  are  similar  to 
those  on  Fe(llO).  The  MEPs  on  those  two  surfaces  are  also 
found  to  be  similar.  Eigure  8  shows  the  first  dehydrogenation 
step  on  Ee3Si(110):  H2S  rotates  in  its  molecular  plane,  bends 
one  H— S  bond  toward  the  surface,  and  then  breaks  off  one  H 
atom  on  a  nearby  Fe  atom.  This  mechanism  is  exactly  the  same 
as  that  on  Fe(l  10).'^  The  second  dehydrogenation  step  to  form 
S  from  SH  also  involves  the  same  bond-breaking  configuration 
for  both  pure  Fe  and  Fe3Si(110).  If  we  look  at  the  Fe3Si(110) 
surface,  we  see  that  the  Fe  atoms  on  Fe3Si(110)  form  a  local 
environment  that  resembles  Fe(l  10).  The  ensemble  size  required 
for  reaction^'  here  seems  to  be  only  4  atoms  of  Fe.  This  similar 
environment  leads  to  similar  barriers  and  pathways  for  H2S 
dehydrogenation . 


-0.8-1  ■ - 1 - . - 1 - . - 1 - . - 1 

0  12  3 

Reaction  Coordinate  (A) 


Figure  7.  PAW-DFT-GGA  minimum  energy  path  for  CO  dissociation 
on  FeAl(lOO)  and  the  stmctures  for  the  initial,  transition  state,  and  final 
states.  Fe  atoms  are  in  dark  gray,  A1  in  light  gray,  O  in  white,  and  C 
in  black. 
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Figure  8.  PAW-DFT-GGA  minimum  energy  path  for  H2S  dissociation  on  FejSilllO)  and  several  intermediate  structures  along  the  path  (the  upper 
images  are  side  views,  and  bottom  are  top  views).  Fe  atoms  are  in  dark  gray,  Si  in  light  gray,  S  in  black,  and  H  in  white. 
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Figure  9.  PAW-DFT-GGA  minimum  energy  path  for  FI2S  dissociation  on  FesSiClOO)  and  several  intermediate  structures  along  the  path.  Fe  atoms 
are  in  dark  gray.  Si  in  light  gray,  S  in  black,  and  FI  in  white. 


The  surface  of  Fe3Si(100)  is  comprised  of  half  Fe  and  half 
Si  atoms,  and  their  lateral  sub-lattices  interleave.  This  affects 
the  Ft2S  dehydrogenation  mechanism.  Figure  9  shows  the  MEP 
for  the  first  dehydrogenation  step.  We  see  that  Ft2S  adsorbs 
parallel  to  the  surface  and  on  top  of  an  Fe  atom  in  the  initial 
state.  Then  the  H2S  molecule  moves  toward  a  nearby  Si  atom 
and  points  the  Ft  atom  of  one  H— S  bond  to  the  Si  atom  at  the 
transition  state.  The  Ft  atom  then  breaks  off,  moves  on  top  of 
the  Si  atom,  and  stays  there.  The  HS  group  moves  back,  passing 
on  top  of  the  Fe  atom  and  then  diffusing  to  the  bridge  site. 
This  pathway  is  Si-mediated  and  thus  very  different  from  that 
on  Fe(lOO),  even  though  the  barrier  is  not  significantly  different. 
After  the  first  step,  HS  adsorbs  asymmetrically  at  the  bridge 
site,  resulting  in  the  H  atom  of  the  HS  group  being  very  close 
to  a  nearby  Fe  atom  (Figure  9).  So  it  is  understandable  that  the 
second  dehydrogenation  step  is  found  to  be  mediated  by  that 


Fe  atom.  Comparing  Fe3Si  with  pure  Fe  for  H2S  dissociation, 
we  see  that  Fe3Si  is  unfortunately  as  susceptible  as  Fe  to  initial 
H2S  attack. 

FeAl(lOO)  is  terminated  with  Al,  and  FeAl(l  10)  does  not  have 
a  local  surface  environment  (reaction  ensemble  size)  similar  to 
Fe(llO).  So  we  expect  different  pathways  for  H2S  dissociation 
on  FeAl  surfaces.  Figure  10  shows  the  MEPs  for  the  first 
dehydrogenation  step  on  FeAl(llO)  and  on  PeAl(lOO).  Even 
though  H2S  starts  with  a  similar  geometry  on  EeAl(l  10)  as  on 
Fe(llO)  and  Fe3Si(110),  the  dissociation  pathway  is  distinctly 
different.  Unlike  on  Ee(llO)  or  Ee3Si(110),  where  H2S  rotates 
in  its  molecular  plane,  instead  H2S  bends  out  of  the  molecular 
plane  toward  the  surface.  However,  as  on  Ee(llO)  or  Fe3Si- 
(110),  where  H2S  breaks  off  the  first  H  atom  on  top  of  an  Fe 
atom,  H2S  also  uses  an  Fe  atom  instead  of  Al  to  dissociate  the 
H— S  bond  on  FeAl(llO).  This  is  also  the  case  for  FeAl(lOO), 
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Figure  10.  PAW-DFT-GGA  minimum  energy  paths  and  several  intermediate  structures  along  the  paths  for  H2S  dissociation  on  FeAl(llO)  (a)  and 
FeAl(lOO)  (b).  Fe  atoms  are  in  dark  gray,  A1  in  light  gray,  S  in  black,  and  FI  in  white. 


where  Fe  is  in  the  subsurface.  In  Figure  10b,  we  see  that  H2S 
bends  one  of  the  H— S  bonds  toward  a  subsurface  Fe  atom,  and 
the  FI— S  bond  breaks  by  transferring  the  H  atom  to  the 
subsurface  Fe  atom.  The  second  dehydrogenation  step  on  FeAl- 
(100)  follows  the  same  pathway  as  the  first  one,  only  with  a 
barrier  that  is  ~0.4  eV  higher.  This  indicates  that  the  HS  group 
may  be  a  stable  intermediate  on  FeAl(lOO)  that  could  be 
observed  experimentally. 

3.4.  Discussion  of  the  CO  and  H2S  Surface  Chemistry  on 
Fe  and  Alloy  Surfaces.  From  Table  4,  we  see  that  H2S  is  a 
more  aggressive  attacker  of  iron,  as  shown  by  the  very  small 
dissociation  barriers  of  H2S  on  Fe.  Interestingly,  the  closest- 
packed  Fe(llO)  surface  turns  out  to  be  more  efficient  than  the 
more  open  Fe(lOO)  surface  in  the  dehydrogenation  of  H2S.  Since 
the  initial  deposition  of  sulfur  is  kinetically  facile,  the  sulfidation 
process  is  probably  controlled  by  other  factors.  Sulfidation  of 
Fe  films  shows  that  the  kinetics  follows  a  parabolic  curve, 
indicating  the  reaction  is  solid-state  diffusion  controlled  (by  Fe 
ions).^^’^^  Accordingly,  the  bulk  diffusion  barrier  is  much  higher 
than  the  H2S  dissociation  barrier.  Grabke  and  co-workers  have 
observed  that  a  small  concentration  of  H2S  in  a  CO/H2  gas 
mixture  can  retard  the  metal  dusting  process. This  is  likely 
to  be  due  to  the  facile  deposition  of  S  atoms  that  block  the 


surface  sites  for  CO  adsorption  and/or  dissociation  (i.e.,  S 
poisons  the  surface).  Our  prediction  of  strong  bonding  between 
Fe  and  S  (~6  eV)  is  consistent  with  this  idea. 

CO  dissociation  on  Fe(lOO)  is  found  to  be  ~0.4  eV  lower 
than  on  the  closest-packed  Fe(l  10)  surface,  indicating  that  more 
open  surfaces  will  lead  to  smaller  barriers  for  CO  dissociation. 

50  we  expect  that  surface  defects  and  steps  on  actual  steel 
surfaces  will  enhance  CO  dissociation.  This  agrees  with  the 
observation  that  after  argon  ion  sputtering,  Fe(llO)  becomes 
more  efficient  in  breaking  the  CO  bond.^® 

The  surfaces  of  the  intermetallic  compounds  FeAl  and  Fe3Si 
provide  a  way  to  investigate  the  effect  of  alloying  on  the  Fe 
surface  chemistry  of  CO  and  H2S.  Again,  our  goal  was  to 
explore  if  alloying  could  raise  the  barriers  to  molecular 
dissociation  such  that  initial  attack  of  steel  surfaces  could  be 
inhibited  through  forming  surface  alloys.  From  Table  4,  we  see 
that  use  of  FeAl  surfaces  increases  the  dehydrogenation  barriers 
of  H2S  up  to  0.4  eV.  Unfortunately,  the  CO  dissociation  barrier 
is  smaller  on  FeAl  surfaces,  especially  on  FeAl(lOO)  [0.7  eV 
lower  than  on  Fe(lOO)].  On  FesSi  surfaces,  the  trend  is  reversed. 
The  resistance  against  initial  CO  attack  is  increased,  but  the 
small  H2S  dissociation  barriers  are  basically  unchanged  on  Fes- 

51  surfaces  compared  with  Fe  surfaces.  Our  work  indicates  that 
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A1  is  not  a  good  choice  for  an  alloying  element  in  terms  of 
preventing  initial  CO  attack  and  that  H2S  resistance  is  difficult 
to  achieve  by  alloying.  If  carburization  is  the  only  concern,  then 
adding  Si  to  the  steel  surface  could  be  helpful. 

4.  Summary  and  Conclusions 

Using  periodic  density  functional  theory  within  the  general¬ 
ized-gradient  approximation  to  electron  exchange  and  correla¬ 
tion,  we  have  studied  CO  and  H2S  adsorption  and  dissociation 
on  low-index  Fe  alloy  surfaces.  We  have  investigated  chemistry 
on  the  (110)  and  (100)  surfaces  of  Fe,  FeAl,  and  Fe3Si.  We 
find  that  CO  dissociates  via  a  lying-down,  bridging  transition 
state  on  all  (110)  and  (100)  surfaces  examined.  FI2S  dissociates 
via  an  on-top-of-Fe  path  on  (110)  surfaces.  The  similar  local 
reaction  ensemble  size  of  Fe  atoms  on  FesSill  10)  compared  to 
Fe(llO)  leads  to  analogous  pathways  and  barriers  for  H2S 
dissociation.  Interestingly,  subsurface  Fe  participates  in  CO  and 
H2S  dissociation  on  FeAl(lOO).  We  found  a  precursor  state  of 
CO  on  FeAl(lOO)  with  an  extremely  low  C— O  stretching 
frequency  and  a  low  dissociation  barrier.  FeAl  is  predicted  to 
be  more  resistant  to  initial  H2S  attack,  but  more  vulnerable  to 
initial  CO  attack  than  Fe.  FesSi  is  predicted  to  be  more  resistant 
to  initial  CO  attack  than  Fe,  but  as  susceptible  as  Fe  to  initial 
H2S  attack.  Neither  Si  nor  A1  alloying  with  Fe  can  simulta¬ 
neously  protect  steel  against  the  initial  stages  of  chemical 
degradation  by  CO  and  H2S.  However,  future  investigations  of 
the  subsequent  solid-state  chemistry  (diffusion  and  compound 
formation)  to  see  if  A1  or  Si  could  inhibit  the  subsequent 
degradation  of  steel  would  be  worthwhile.  If  this  next  avenue 
also  proves  ineffective,  an  alternative  approach  to  protecting 
steel  is  to  develop  coatings  that  can  form  protective  oxide  scales 
and  that  do  not  fracture  easily.  We  are  pursuing  such  ideas  at 
present. 
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